Effective actions for spin ladders. 
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We derive a path-integral expression for the effective action in the continuum limit of an AFM 
Heisenberg spin ladder with an arbitrary number of legs. The map is onto an 0(3) nonlinear cr- 
model (NLcrM) with the addition of a topological term that is effective only for odd-legged ladders 
and half-odd integer spins. We derive the parameters of the effective NLcrM and the behaviour of 
the spin gap for the case of even-legged ladders. 
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Triggered by the discovery of high-T c superconduc- 
tivity and by the fact that, in the strong-coupling 
regime, the Hubbard model maps onto an AFM Heisen- 
berg model quantum spin systems have become in 
recent years of great theoretical and experimental inter- 
est H . Independently from the possible connections with 
high-T c materials, some years ago Haldane j| put for- 
ward the conjecture that the value of the spin should 
discriminate dramatically between half-odd integer and 
integer one-dimensional spin chains: the former should be 
gapless, with power-law decay of spin-spin correlations, 
while the latter should be characterized by exponential 
decay of correlations and hence by a nonvanishing spin 
gap. It turns out that, in the continuum limit, an 
AFM Heisenberg chain is described by an 0(3) nonlin- 
ear cr-model (NLcrM) B with the addition of a topolog- 
ical term (a Pontrjagin index multiplied by a coef- 
ficient 9 — 2irS, where S is the value of the spin. As 
shown, e.g. by Shankar and Read 0, it is precisely this 
term that makes the half-odd integer spin chain (where 
9 = 7r mod 2ir) massless, while it is ineffective for integer 
spin chains (where 9 = mod 2tt), which are gapped. 

Extensions of these results to D=2 turned out to be 
rather disappointing, to the extent that it has been 
proved in a convincing way [8-12] that, at least for smooth 
configurations (see however |1 for the case of a singular 
"hedge-hog" field configuration), the topological term is 
absent irrespective of the value of the spin and also of 
the topology of the lattice. 

Quite recently, however, it came as a surprise when 
theoretical (mostly numerical) and subsequently experi- 
mental results showed that spin ladders, obtained by an- 
tiferromagnctically coupling a finite number of chains, 
show drastically different behaviours that are basically 
determined by the value of the spin and by the number 
of legs in the ladder. Coupled spin chains had already 
been studied previously 113] mainly as models of single 



chains with spin higher than the minimum value S — 1/2. 
A systematic study of the appearance of a spin gap in 
even-legged ladders and of the possible onset of super- 
conductivity upon doping with holes was initiated later 
on by Dagotto et al. [B who made for the first time 
quantitative predictions on both the spin gap and on the 
pair correlation function. Odd-legged ladders are instead 
gapless. This "even-odd" conjecture was formulated for 
the first time by Rice et al. ]l5| , and is very much remi- 
niscent of the already mentioned Haldane's "integer-half 
integer" conjecture Q for single spin chains. According 
to this "even-odd" conjecture, undoped integer-spin lad- 
ders are gapped while half-odd integer spin ladders are 
gapped when the number of legs is even, gapless when it 
is odd. These results have by now strong support, both 
theoreticall y [l6| and experimentally (l7j (for a recent 
review see |18|). The situation seems to be not so sim- 
ple for conducting (i.e. hole-doped) ladders. A different 
behaviour between even and odd-legged ladders is pre- 
dicted only for very strong intrachain repulsions ]fl9|| . 
On the contrary, in the weak coupling regime, analytical 
fl9| , p0f and numerical Q studies suggest that d-wave 
interchain pairing correlations become dominant over AF 
fluctuations in both the two and the three-legged ladder, 
due to the simultaneous presence of gapless and gapped 
spin modes. So, the situation is richer than in the case of 
single conducting chains, that display a Luttinger liquid 
type (THJH] behaviour. 

Following mainly the ideas of Haldane 0] and Af- 
fleck H on the role of the topological term in spin chains, 
considerable theoretical effort has been recently devoted 
to the investigation of the existence and the role of such a 
term when a finite number of chains are coupled to form 
a ladder (the possible topological origin for the different 
behaviour of even and odd leg ladders has been suggested 
for the first time in (23|). In particular Senechal has 
given a derivation of the NLcrM continuum limit of a 
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two- leg ladder using a coherent-state path-integral |2J] 
expression for the partition function, while Sierra [26] 
has employed an (operator) Hamiltonian approach fol- 
lowing closely Affleck's mapping of the Heisenberg chain 
onto the NLcrM. 

What we present here is a path-integral approach that 
is however different from Senechal's. It allows us to re- 
produce Senechal's results for the two-leg ladder, but also 
to generalize them to ladders with an arbitrary number 
of legs and to obtain a clear understanding of why the 
topological term is absent for even- legged ladders. 

The Hamiltonian for a ladder system with n\ legs of 
lenght N is defined by: 



n, N 



a—1 i—1 
J'a,a+lSa{i) ' S a+1 (i)} . 



(1) 



The only condition we shall impose on the coupling 
costant J a and J' a ii is that the classical minimum of 
the Hamiltonianian (Fy) be antiferromagnetically ordered. 
The partition function for the Hamiltonian (jl]) in a path- 
integral representation which makes use of spin coherent 
states J2f| is given by: 

Z(/3) = J [£>«] exp J is w[«a(*, 7")] - f dTH(r) 

(2) 

where u>[Sl a (i, r)] is the Berry phase factor coming from 
the exponentiation of the overlap between coherent states 
at nearby timeslices f27|| , while H(t) is obtained by re- 
placing in the Hamiltonian (|l|) the operator S a {i) by the 
classical variable stt a (i,T). Explicitly: 



u[CI(t)] = dr0(r)(l- cos 0(r)) (3) 
Jo 

h(t) = {Jas 2 n a {i, t) ■ n a {i + 1 , t) 

a,i 

+ j' aM+lS 2 n a (t,T)-n a+1 (i,T)} . (4) 

The Berry phase measures the area enclosed by the 
path fl(r) = (sin#(Y) cos </>(t), sin 9(t) sin 4>(t), cos 9{t)) 
on the unit sphere. We shall proceed now along the lines 
of Ref. ||. We assume that short-range antiferromag- 
netic correlations survive at the quantum level, so that we 
consider the dominant contribution to the path-integral 
as coming from paths described by: 



tl a {l 1 T) = {-l) a+l ^T) 1 



1/2 



l a {i,T) 



The fluctuation field l a (i) is supposed to be small: 
\l a (i)/s\ « I and the field <p(i,r) slowly varying. We 
are then allowed to make an expansion up to quadratic 

order in I, <fi and 0. The parametrization of the classical 
fields in (phis different from that employed in Senechal's 
approach |24| to the two-leg ladder, and allows for gener- 
alization of the path- integral approach to ladders with an 
arbitrary number of legs. In particular, the fact that the 
field 4> is chosen to depend on the site index i along the 
legs, but not on the index a labeling the sites along the 
rung, reflects the assumption that the correlation length 
be much greater than the total width of the ladder, i.e. 
£ >> nia, where £ is the staggered spin-spin correla- 
tion length and a the lattice spacing. This is well sup- 
ported by several numerical works (see for instance |2§| l 

and |^]). The constraint & a {i) = 1 implies cj) (i) = 1 
and <j)(i) ■ l a (i) = 0. As far as the intra-leg term of H{t) 
is concerned, everything proceeds as for the continuum 
limit of the I-dimensional Heisenberg chain (see for 
a detailed calculation) and hence we write just the final 
result: 

J a s 2 n a (i,T) ■n a (i + l,r) M 

dx ^ ( * 2 ^° Ja) i\x, r) + 2 Y J a \l a (x, r)] 2 j . (6) 
For the inter-leg term we have instead: 



Sl a (i,r) ■ O a+ i(i,r) w 

\l a (i,T)\ 2 \l a+ i(i,r)\ 2 l a (i,r) -l a+1 (i,T) 
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The term H(t) in the action has then the continuum 
limit: 



H(r) 
I 



dx 



(s 2 Y Ja)4> (X, T) + Y la(x, T)L a ^l b {x, T) 



a . b 



where L a ^ is the same matrix defined in Ref. [|26|| , i.e.: 



Jl 

J a,b 



(J'a,a-1 = J 'a-i.a an d J[ fl = J' nu n l+1 = in the formula 
above). We have finally to evaluate the Berry phase term. 
We now need the formula for the variation of the Berry 
phase lo[4>] upon a small change 8<j> [E7fl: 



5lo 



drdcj) -(0x0) 



(9) 
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We have therefore, at leading order: 
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i,a i,a 

22 dr (cp(i, t) x cp(i, t)) ■ l a (i, t) 



r&] + J2 dT(j>(i,T)x i>(i,T))-i a (i, T ) . (10) 



The first term in the last equation is the topological term. 
We shall return to it in a while, after having integrated 
out the fluctuation field l a (i,r): 



D[l] exp J dr J dx 



~n l a( X > T)L a , b l b (x, t) 
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dr / dx 
o J 2 



(11) 



The integration over the field l a (i) gives therefore the 
kinetic term of the NLerM. For the topological term: 



r[0] = s y>iry>i)v#(*,T)], 



(12) 



we can simply observe that, since the field <p(i,r) does 
not depend on a, one has just the same topological term 
as for the chain times a factor which is zero for even rii 
and 1 for odd n;: 

T[<f>] = ( -L$n dT $ dx <l> ■ (0 x <t>) m odd ( 13 ) 
I ni even 

(6 = 2irs in the last equation). In the case of the 2D (infi- 
nite) lattice, cancellation of the topological term results, 
as is well known [8-12], from taking the continuum limit 
also along the direction of the rungs, and holds at least 
for smooth field configurations. Here it is instead a direct 
consequence of the assumption that £ >> n/a, which lies 
at the heart of the parametrization chosen in Eq. (|^). 
The staggering of the field (f> along the rungs leads then 
to the result that the topological term survives only for 
odd-legged ladders, being of course significant only for 
half-odd integer spins. Putting everything together we 
found that the antiferromagnetic Heisenberg ladder sys- 
tem is mapped onto a (1+1) NLerM with the euclidean 
Lagrangian: 



C = < 



M ±^+vJ 2 ) +£ 0-(0X0) n z odd 



2y 



ni even 



(14) 



where the NLerM parameters are defined by: 

1/2 
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Ylb,C ^b,C , 



1/2 



(15) 
(16) 



Let us remark that the NLerM velocity we have obtained 
coincides with the spin- wave velocity (see J2(| for the cal- 
culation of the spin-wave velocity in our model Hamilto- 
nian). In the particular case n; = 2 we obtain: 



-(1 + J'/2J) 1 / 2 



v s = 2sJ{\ + J'/2J) 1/2 



9 



(17) 
(18) 



These values of the parameters concide with the ones 
obtained by Senechal in H] for the two- leg ladder. 
They are different, however, from the parameters ob- 
tained by Sierra through his Hamiltonian mapping to the 
NLerM He has for the n ; -leg ladder: 



v s = s 



2 E JaL b,l ~ \ 5n i 



o J2 a Jg _ r 



1/2 



( 2 J2b,c L b,l 



1/2 



(19) 



(20) 



where 5 ni is equal to 1 for odd ni and for even n;. 
For even ni his g is smaller, while v s is greater, by a 
factor of y/2 as compared with ours. The two results 
coincide (and coincide with Afflecks's as well, as they 
should) in the limit of a single chain, i.e for ni = 1. 
They are again different for odd ni > 1. We argue that 
this is due to the introduction, in Ref . j[6) , of additional 
massive fields. The latter are decoupled from the er-modcl 
field and are ultimately neglected in the effective theory 
developed in |2(|. Nonetheless, they appear to affect the 
actual values of the parameters of the effective er-model. 
The agreement of our approach and of that of Ref. |2lf | 
for ni — 1 (when there is no room for massive fields) 
seems to give support to our conjecture. 

In any case the general picture for the behaviour of the 
spin-ladder is the same as in [p6[ . The ladders with an 
even number of legs are mapped to a (1+1) NLerM with- 
out topological term and are therefore gapped, while the 
ones with an odd number of legs are described by a (1+1) 
NLerM with the topological term, and are gapless for 
half-odd integer spin. In the case J' a a+l — J', J a = J 
and J' << J we have g ~ 2/sn ; ; the coupling costant 
g gets therefore smaller and smaller when the number 
of legs increases and the NLerM enters in a weak cou- 
pling regime. In this regime we may use the formula 
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A ~ exp(— 2n/g) |3(j] to estimate the spin gap for even- 
leg ladders. We obtain then A ~ exp(— nsni); the gap 
decreases therefore with the number of legs, as observed 
in numerical simulations p8j , p9| and as expected from the 
fact that in the limit ni — ► oo the difference between lad- 
ders with even and odd number of legs must disappear. 
Let us also remark that the ratio £,/nia ~ e^jp(nsni)/ni; 
the condition £ << uja we supposed at the beginning 
of our calculation seems therefore to be satisfied (in a 
self-consistent way) better and better when the number 
of legs increases (at least for J' « J). In the opposite 
regime J' » J the coupling costant becomes strong 
since g ~ (J'/J) 1 ^ 2 and we can therefore estimate the 
spin gap as A = v s g in the strong coupling regime 0. 
For two legs we get A ~ J' when J' >> J, a result which 
also agrees with what found in literature through other 
techniques |l8[. 



Remark. Just after completion of this work a preprint 
by G. Sierra has appeared pl[ , in which the mapping of 
spin ladders to the NLcrM is reviewed also within a path- 
integral formalism. The values for the NLcrM parameters 
found by Sierra agree with ours. 
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